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ABSTRACT 

We discuss wide-field near-IR imaging of the NGC1333, L1448, L1455 and Bl 
star forming regions in Perseus. The observations have been extracted from a much 
larger narrow-band imaging survey of the Taurus- Auriga-Perseus complex. These H 2 
2.122/im observations are complemented by broad-band K imaging, mid-IR imaging 
and photometry from the Spitzer Space Telescope, and published submillimetre CO 
J=3-2 maps of high-velocity molecular outflows. We detect and label 85 H2 features 
and associate these with 26 molecular outflows. Three are parsec-scale flows, with a 
mean flow lobe length exceeding 11.5'. 37 (44%) of the detected H2 features are as- 
sociated with a known Herbig-Haro object, while 72 (46%) of catalogued HH objects 
are detected in H 2 emission. Embedded Spitzer sources are identified for all but two 
of the 26 molecular outflows. These candidate outflow sources all have high near-to- 
mid-IR spectral indices (mean value of a ~ 1.4) as well as red IRAC 3.6 /itm-4.5 /mi 
and IRAC/MIPS 4.5 ^m- 24.0 ^m colours: 80% have [3.6]-[4.5]>1.0 and [4.5]-[24]>1.5. 
These criteria - high a and red [4.5]-[24] and [3.6]-[4.5] colours - are powerful discrim- 
inants when searching for molecular outflow sources. However, we find no correlation 
between a and flow length or opening angle, and the outflows appear randomly orien- 
tated in each region. The more massive clouds are associated with a greater number 
of outflows, which suggests that the star formation efficiency is roughly the same in 
each region. 

Key words: stars: formation - infrared: stars - ISM: jets and outflows - ISM: 
kinematics and dynamics - ISM: individual: Perseus 



1 INTRODUCTION 

The U.K. Infrared Deep Sky Survey (UKIDSS) is a legacy 
survey bei ng conducted at the United Kingdom Infrared 
Telescope l|Lawrence et al.ll2007D . UKIDSS comprises five 
surveys, ranging from very deep, narrow studies of extra- 
galactic fields, to shallow, broad surveys at both low and 
high galactic latitudes. One such survey, the Galactic Plane 
Survey (GPS), will by 2012 observe 1868 deg 2 of the north- 
ern and equatorial galactic plane between -5° < b < +5°, 
in J, H and K down to limiting magni tudes of 20.0, 19.1 
and 19.0, respectively (|Lucas et al.ll2008h . The GPS also in- 
cludes a component covering the Taurus- Auriga-Perseus (T- 
A-P) complex of star-forming clouds, through narrow-band 
H 2 1-0S(1) and broad-band K filters. To date, 140 deg 2 of 
the T-A-P have been observed as part of the GPS. Obser- 
vations of an additional 40 deg 2 are currently being sought 
to complete the proposed survey area, which is shown in 



Fig. [T] The existin g data are availab le from the WFCAM 
science archive (see Lucas et al J 1200a. for details of how to 



access these data) and, as mosaicked images, from a dedi- 
cated database hosted by the Joint Astronomy Centr^B using 
an easy-to-navigate clickable map. The use of the latter is 
recommended for those wishing to quickly access imaging 
data to search for H2 emission features. 

The goal of this paper is to present demonstration sci- 
ence from the T-A-P H2 survey. We therefore focus on just 
the western end of Perseus, comparing the narrow-band 
imaging data with photometry from the Spitzer Space Tele- 
scope legacy programme From Mol ecular Cores to P lanet- 
Forming Disks, also known as c2d (|Evans et al.ll20o"3 ). An 
area covering almost 4 deg 2 (representing five WFCAM 
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Figure 1. An overview of the Taurus- Auriga- Perseus complex. 
The grey-scale shows integrated CO J=l-0 emission from a survey 
of the Milky Way obtained with 1.2 m telescopes in Cambridge, 
MA and on Cerro Tololo, Chile {5 ame. Hartmann fc Thaddeusl 
1200 it) . The contours represent visual extinction, A v , con- 
structed from optica l star-counts using Digitised Sky Survey data 
(Dobashi ct al. 2005). The black lines outline the complete region 
to be observed in H2 2.122 fim emission and broad-band K. 



Figure 2. Colour-magnitude diagram of c2d Perseus data. Black 
dots represent the entire c2d Perseus catalogue; red circles mark 
those identified as YSOs, while blue crosses mark those labelled 
"red" in the catalogue. 



2 OBSERVATIONS 



2.0.1 WFCAM data 



tiles) has been searched for H2 2.122 /im line emission fea- 
tures. The mid-IR data are then used to identify candi- 
date protostars for the molecular outflows found in H2. We 
also search for H2 counterp arts of the known Herbig-Haro 
(HH) o bjects catalogued bv lWalawender. Bally fc Reipurthl 
1 2005 j). and compa r e our i mages with the CO outflow maps 
of iKnee fc Sandelll l|200rj| ) and lHatchell. Fuller fc Richer! 
(2007a|). Our aim is to establish the relative effectiveness 
of optical, near-IR and CO surveys of outflows. 

The Perseus cloud complex spa ns more than 30 pc at 
a distance o f between 200-350 pc (|de Zeeuw et al.l 1 19991 ; 
Cernisl |l990T) The region was first mapped in CO by 



Sargent (1979); molecular line and ex ti nction maps have 



since been presented by I Ridge et al.l l|2006l ). The early 
search for Young Stellar Objects (YSOs) an d protostars 
from IRAS data bv lLadd, Lada fc Myers! d 19931) has recently 
been su perseded by the Sp itzer s urveys of |j0rgensen et all 
l|2006l) . lLada et all \ 2006) and iGutermuth et al. I (|200Sl ). 
Perseus harbours two particularly active star forming re- 
gions, NGC1333 - centre d on the bright and very well stud- 
ied HH complex HH 7-11 dHodapp fc Laddll995l ; iBallv et all 



1 19961 ; IChrvsostomou et alfcOOOh- and IC 348, a r ich, embed- 
ded stellar cluster ( Luhman et al ] Il998l . 120051 ; lLada et al.l 
2006). Here we focus only on NGC1333 and surrounding 
clouds - a region we refer to as "Perseus- West" - although 
the full extent of Perseus will be imaged in H2 as part of the 
T-A-P survey (F ig.D, and in JHK by the UKIDSS Galactic 
Clusters Survey (|Lawrence et alj|2007l ). 



The near-IR wide-field camera WFCAM (|Casali et al.ll2007h 
contains four Rockwell Hawaii-II (HgCdTe 2048x2048) ar- 
rays spaced by 94% in the focal plane. The pixel scale mea- 
sures 0.40". To observe a contiguous square field on the sky 
covering 0.75 deg 2 , a region known as a WFCAM "tile", ob- 
servations at four positions are required. At each of these 
four positions eight exposures are obtained in each filter, 
with a 20 sec exposure time in H2 and a 5 sec exposure 
time in K. The images are shifted by an integral multiple of 
0.2" so that in the final, mosaicked images the pixel scale 
is 0.2"; in this way the sub-arcsecond seeing typically ob- 
tained at Mauna Kea is fully sampled and bad pixels may 
be compensated for. 

All WFCAM data are initially reduced by the Cam- 
bridge Astronomical Survey Unit (CASU) and are dis- 
tributed through a dedicated archive hosted by the Wide 
Field Astronomy Unit (WFAU). The CAS U redu ction 
steps are describ ed in detail by iDve et all (l2006h and 
llrwin et al.l l|2008l ); astrometr ic and photomet r ic calibrations 
are a chieved using 2MASS (|Dve et al.l 120061 ; iHewett et all 
2006). Data are then downloaded in bulk from WFAU 
and combined into tiles using astrometric information 
stored in the file headers. As part of this process, resid- 
ual sky/background structure is removed b y fitting a 
coars e surface to each image (as described by iDavis et al.l 
120071 ). Currently data for 140 deg 2 are available; H2 and 
K-band mosaic tiles are available to all interested par- 
ties from a database at the Joint Astronomy Centre 
(http://jach.hawaii.edu/UKIRT/TAP) as Rice-compressed 
fits files. 
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2.0.2 Spitzer data 

The Spitzer Space Telescope observations discussed in this 
paper were obtained w ith the IRAC and MI PS cameras, 
which are described by iFazio et al.l l|2004 ) and iRieke et al.l 
l|2004 l. respectively. IRAC data in bands 1,2,3 and 4 (at 3.6, 
4.5, 5.8 and 8.0 (im) and MIPS band 1 data (at 24 fj,m) were 
obtained in Sept ember 2004 as par t of the Spitzer legacy 
programme c2d dEvans et all 120031) . The data have been 
presented by iJorgensen et al.l (|2006l ) ; gu aranteed time ob- 
servat ions of NGC1333 are presented by iGutermuth et ail 
(2008). Additio nal data on L1448 were obtained on Febru- 
ary 25, 2005 bv lTobin et ail (|2007l ). 

The mid-IR photometry used here were 
obtained directly from the c2d archive 
(http://peggysue.as.utexas.edu/SIRTF/), which includes 
data extracted from the third delivery of IRAC and MIPS 
data. The procedures used du ring source extrac t ion ar e de- 
scribed by IE vans et~ail (|2005l ) and lHarvev et all (|2007f ). For 
Perseus, the full catalogue of sources detected in at least one 
IRAC band contains over 120,000 objects. IRAC and MIPS 
band 1 (24 /im) data have been combined with 2MASS 
JHK photometry by the c2d group to distinguish candidate 
Young Stellar Objects (YSOs) and proto stars from field 
stars and background galaxies (see also Ijprgensen et al.l 
2006). A spectral index, a = dlog(XF(\))/dlog\, calculated 
between 2.2 /im and 24 /jm, is included in the catalogue 
for all source with at least two photometry points. In this 
paper we over-plot the positions of YSOs and protostars 
identified in this way, as well as the "red", "redl" and 
"red2" sources (targets with a high MIPS band 1 to IRAC 
flux ratio, or those with rising spectral energy distributions 
[SEDs] detected in only MIPS1 or MIPS1 and IRAC band 
4) and "PAH-em" and "star+dust" sources (targets with 
IRAC colours expected for PAH emission features; these 
may include background galaxies but also Class II sources). 
These latter sources are arguably also candidate protostars 
and/or outflow sources. 

Note that both this paper and the c2d catalogue refer 
to all young stars as YSOs, regardless of whether they are 
Class 0/1 protostars or Class II /III young stellar objects. 

Sample SEDs, plotted betwee n 1.2 /im and 24 ^m , are 
presented for all source types by lEvans et al.l l|2005l ). In 
Fig. [2] we plot a Spitzer colour-magnitude diagram for 
sources in Perseus West from the c2d catalogue. There are 
66,535 sources in the catalogue, of which 4,352 are "red" 
sources and 610 YSOs. Clearly, most of the "red" sources, 
though red in the two IRAC bands, have a low 8.0 /im flux. 
In the most deeply embedded objects this may be caused by 
silicate absorption or the intrinsic shape of the SED (Class 
sources have been shown to possess a [3.6]-[4.5]/[5.8]-[8.0] 
ratio that is >1; J0rgensen et al. 2006). Even so, the "red" 
sources may also include a modest population of b ackground 
galaxies l|Harvev et al.|[2006l ; Ijprgensen et al.ll2006h . 



2.0.3 J CM T data 

The H2 emission in outflows is known to be produced in 
molecular shocks formed at the interface between a jet 
and the ambient medium. Thus, the H2 knots and bow 
shocks imaged with WFCAM are expected to be closely re- 
lated to the high- velocity molecular flow lobes often mapped 



at (sub)millimetre wavelengths in rotational CO emission. 
In Figs[3] to [7] we therefore present CO J=3-2 observa- 
tions of high-velocity molecular outflows in Perseus- West 
taken from the literature. These data were acquired at the 
James Clerk Maxwell Telescope (JCMT) and have been pub- 
lished and discussed in de t ail bv|Knee fc Sandelll (|2000l ) and 
lHatchell. Fuller fc Richer! l|2007al ). The JCMT beam width 
measures 14" at the 345.796 GHz frequency of the J=3^2 
transition. Contour plots are presented in units of antenna 
temperature. 

The maps cover only very limited regions in each star- 
forming cloud, though many of the flows discussed below 
have been observed. More extensive maps are being planned 
as part of the Gould Belt legacy survey at the JCMT 
jWard- Thompson et al.1l2007l ). 



3 RESULTS & DISCUSSION 
3.1 H2 imaging of Perseus- West 

The five WFCAM tiles used in this paper cover the west- 
ern portion of Perseus, which includes NGC1333, L1455, 
L1448, and Barnard 1 (Bl). WFCAM H 2 images of regions 
where outflows were detected are shown in Figs (3HS1 The 
positions of protostars and YSOs identified in the Spitzer 
c2d survey are overplotted in each figure. Circles indicate 
sources that were given the "YSO" designation in the c2d 
catalogue, while triangles mark sou rces classified as "red" , 
"red2" , "PAH-em" , or "star+dust" . |j0rgensen et all (|2007f ) 
have compared the c2d data with SCUBA observations at 
850 /im. They present a catalogue of 49 "embedded" YSOs; 
all are MIPS detections that are either red in [3.6]- [4.5] and 
[8.0]-[24] colours or are located within 15" of SCUBA peaks. 
They also include a few concentrated SCUBA cores with no 
MIPS counterparts. We label these sources as "YSO" with 
the numbers from J0rgensen et al.'s catalogue (their Table 
3) and list those sources located in the Perseus- West region 
in Table □ 

The locations of Herbig-H a ro (HH ) objects found by 
IWalawender. Bally fc Reipurthl (|2005al ) in their wide-field 
optical imaging, as well as HH objects that were previously 
known, are also marked on our figures. Note that these po- 
sitions are often rather vague, representing the location of 
what is usually an extended object or diffuse feature. 

In the online Appendix we present UKIRT- 
WFCAM/Spitzer-IRAC colour-composite images of 
each region, and discuss the outflows marked in Figs. |3HH] 
in some detail, noting in particular the association of H2 
emission-line features with HH objects, molecular CO out- 
flows, dense cloud cores and the embedded source identified 
through Spitzer photometry. H2 knot parameters (flow 
position angles, flow lengths, associated HH objects, etc.) 
are listed in Table IA11 where we also give the embedded 
YSO that is likely to be driving each flow. Many of these 
knots or groups of knots are known HH objects or parts of 
outflows imaged in less extensive surveys. We therefore also 
list the names of associated HH objects and previous H2 
outflow names in Table IA1I 
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Figure 3. Logarithmically-scaled H2 2.122/jm image of Bl with the locations of HH objects marked with crosses, Spitzer-identified 
YSOs marked with circles and Spitzer-ide ntified red sources marked with triangles. j0rgensen embedded YSOs are labelled. Contours 
of high-velocity CO J=3-2 emission (from lHatchell. Fuller k. Richerll2007ari are over-plotted: integrated antenna temperature is plotted 
from 9 to 15 km s — 1 (dashed lines) and -5 to 4km s — 1 (full lines). Contour levels are 1,2,4,8 Kkm s" 1 . H 2 features are numbered for 
reference (see Table [ATJ . Outflows are indicated with arrows. 
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Figu re 4. H2 2.122 fim imag es of two adjacent regions in the Bl ridge. Contours of high-velocity CO 3-2 emission (from 
lHatchell. Fuller fc Richer 2007a) are over-plotted from 9 to 15km s 1 (full lines) and from -5 to 4km s 1 (dashed lines). Contour 
levels are 0.5,1,2 Kkm s _1 (left panel) and 1,1.5,2 Kkm s — 1 (right panel). See Fig|3]for further details. 



3.2 On the merits of combining optical, infrared, 
and submillimetre observations of outflows 

Overall, the combined WFCAM H2 imaging, Spitzer imag- 
ing and photometry, and CO outflow observations provide 



a clear and rather complete picture of dynamic activity in 
Perseus- West. Our WFCAM observations demonstrate that 
even shallow H2 imaging is an effective tool for finding out- 
flows from the youngest YSOs. 

Of the 26 molecular outflows in our WFCAM tiles, 24 



H 2 survey of Taurus-Auriga-Perseus 5 




Figur e 5. H2 2. 122 fim image of L1455 with L1451 to the west. Contours of high-velocity CO 3-2 emission ffrom lrYatchell, Fuller fc Richer! 
2007a) are over-plotted from 6 to 15 km s 1 (full lines) and from -5 to 3km s 1 (dashed lines). Contour levels are 0.5,1,2 Kkm s _1 . See 
Fig|3]for further details. 



of them (92 %) are driven by source s identified as embed- 
ded YSOs bv ljorgensen et al.l \200H ): only 9% of the 85 H 2 
features have no obvious CO flow counterpart or j0rgensen 
YSOs. The two outflow sources not identified by j0rgensen 
et al. as embedded YSO, LkHa327 and L1448-IRS1, were 
probably missed because of saturation in the Spitzer bands. 
Neith er source coincides with a SCUBA core (jHatchell et al.l 
120051 ), so these bright sources have probably already cleared 
their cores. Of the 38 j0rgensen YSOs in ou r field, 24 (63%) 
drive outflows. Clearly, the criteria used by |j0rgensen et all 
<|2007h (MIPS/IRAC colours [3.6]-[4.5]>l and [8.0]-[24]>4.5 
and/or associated SCUBA cores) is not only very good at 
identifying YSOs, but also at identifying YSOs that drive 
molecular outflows. 

Shock-excited outflow features can of course be de- 
tected at both optical and infrared wavelengths, by virtue 
of their forbidden emission lines ([OI], [SII], [OIII], etc.), 
hydrogen recombination lines (Ha), or - in the near-IR 
and mid-IR - molecular hydrogen ro-vibrational lines. How- 
ever, the longer-wavelength transitions will be more effec- 
tive in regions of high extinction. Also, these lines derive 
from very different gas components: the optical lines typ- 
ically trace hot, dense, partially-ionised, high-velocity gas 
(T~10,000K; flow velocities approaching 100-200 km s _1 ) 
while the H2 lines trace low-excitation, shocked molecular 
gas at much lower velocities (T~500-2,000 K; v~10-50 km 
s _1 ). It is theref ore perhaps not surprising that, of the 158 
HH objects from IWalawender. Bally fc Reipurthl (|2005al ) in 
our field, only 72 (46%) of them were detectable in H2 



2.122 /im emission. Similarly, only 37 (44%) of the 85 la- 
belled H2 features are associated with an HH object. Many 
of the HH objects are faint in the optical, so in a few cases 
the lack of a detection in near-IR H2 emission may be due 
to the modest sensitivity of the WFCAM survey. However, 
many of the HH objects are located around the periphery 
of each star-forming cloud (see for example the comparison 
of HH and H2 positions in Fig. [7} , where the molecular gas 
density may be too low to facilitate H2 excitation. Clearly, a 
colour-composite image comprised of optical ([SII] or Ha), 
near-IR 2.12 fim and mid-IR 4.5 fim imaging would be ideal 
for illustrating the high-excitation atomic/ionised (HH), hot 
molecular (ro-vibrational H2), and warm molecular (pure- 
rotational H2) flow components simultaneously, as well as 
the colours of white/blue foreground stars through to the 
reddest protostars. Such a comparison is obviously now pos- 
sible with ground-based optical imaging, WFCAM survey 
data, and Spitzer observations, for a number of well-known 
star forming regions. 



3.3 General H2 flow characteristics 

Of the 26 H2 outflows identified in column 2 in Table (TJ 
20 (77%) exceed an arcminute in mean length; three are 
parsec-scale flows, with a mean lobe length exceeding 11.5' 
(the outflows driven by YSOs 2, 5, and 20c). Longer flows 
might be expected from more evolved sources, and indeed 
YSO 20c and YSO 5 do have flat spectral indices. However, 
YSO 2 has a steeply rising SED (a ~ 2.3). This high value of 
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Figure 6. H2 2.122 (im image of L1448. Contours of high-velocity CO 3-2 emission (from IrTatchell. Fuller fc Richedl2007aT i arc over- 
plotted from 8 to 15km s _1 (full lines) and from -5 to Okm s _1 (dashed lines). Contour levels are 1,2,4,8 Kkm s _1 . Sec Fig|3]for further 
details. 



a may in part be due to excess emission from lines in some of 
the Spitzer bands (discussed below). However, when we plot 
flow length against Spitzer spectral index we find no correla- 
tion between these two parameters. This suggests that even 
the most deeply embedded H2 jet source are old enough to 
produce parsec-scale flows. It seems likely that jet length, 
when derived from H2 observations, is more sensitive to the 
environment (specifically the presence of ambient molecular 
gas that may be shock-excited into emission) than the "age" 
of the jet central engine. 

Can we say anything about flow opening angles from our 
sample of jets in Perseus- West? High- resolution CO maps of 
outflows seem to indicate that Class sources produce more 
collimated flows than the i r more evolved Class I counter- 
parts (e.g. iLee et al.|[2002l ; lArce fc Sargent! [2006). H2 emis- 
sion features (bow shocks) often envelope the outer edges of 
CO flow lobes, so one might expect to see a similar trend 
in H2 data. However, shock-excited H2 emission has a very 
short cooling time (of the order of a few years), so H2 fea- 
tures are often short-lived and confined to a small percentage 
of the flow surface area. Moreover, close to the central source 
H2 emission features tend to be associated with the central, 
on-axis jet rather than the swept-out cavity walls. Defining 
flow opening angles from H2 images is therefore rather dif- 



ficult. Even so, in Table [T] we list the mean opening angle, 
6, for well-defined flows with multiple, bright H2 emission 
features. 6 is measured from cones in each flow lobe that 
are centred on the outflow source and which include all H2 
line-emission features. A plot of 9 against source spectral 
index a (not shown) reveals no correlation between these 
two parameters. At best, we can only say that most of the 
flows in Perseus- West appear to be well collimated, as is 
expected for jets from young, embedded (high-a) sources. 
The opening angles listed have a mean of 12.5° and a stan- 
dard deviation of only 4.1°. The associated sources have a 
relatively broad range in a: -0.87 to 3.79. 

Many of the H2 jets in Perseus appear to be curved or 
possibly even precessing. If both lobes curve in the same di- 
rection (H 2 features 2-4 in the Bl ridge and 76/77 in L1448 
are good examples) the flow curvature may be due to the 
motion of the source through the ambient medium. In such 
cases the source would need a proper motion that is a sizable 
fraction (at least a few percent) of the jet velocity. Source 
tangential velocities of the order of 1-10 km s _1 would be 
required. Alternatively, the apparent curvature may simply 
be due to excitation of H2 along only one side of an oval 
cavity, perhaps because of enhanced ambient gas densities 
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Figure 7. H2 2.122 fim image of NGC1333 . Open squares mark SCUBA cores without c2d-identified YSOs; the other symbols are as in 
Fig|3] High- velocity CO 3-2 contours (from lKnee fc SandelfeOOOh are over-plotted, from 12 to 18 km s 1 (full lines) and from -5 to 3 km 
s _1 (dashed lines). Contour levels are 1.5,3,6,12 Kkm s _1 . 



or an ambient density gradient across the width of the flow 
lobe. 

Precession or jet " meandering" has been identified in 
a number of flows (e.g. lLopez et a l, 1995; Heat hcote et al.l 
ll996l ; lDavis et al.lll997l ; lHodapp~ al. 2005). Precession may 
result from the misalignment of binary orbital and circ um- 
binary disk axes ijTerauem et al. I ll999l : iBate et al.llifJOOh , al- 
though the disk should be forced into alignment with the 
orbital plane in roughly one precession period, which should 
be about 20-tim es longer than the binary orbital period 
|Bate et al.ll2000l ). Given this time limitation, jet precession 
is expected from only the most embedded and therefore po- 
tentially very youngest sources. Notably, YSOs 11, 14 and 
38, the three precessing jet candidates in Perseus- West, all 
have high values of a. 

The H2 flows from YSOs 11, 14 and 38 appear to have 
precessed through roughly half a turn. The dynamical age 



of each flow should thus be a fair estimate of the prece ssion 
period, to within a factor of ~2 l|Hodapp et al ] 120051 ). For 
a canonical jet velo city of 100 km s -1 (con sistent with jet 
proper motions, e.g. lEisloffel fc MundtJ li994). from the flow 
lengths listed in Table [TJ flow ages of 7300 yrs, 1400 yrs 
and 2600 yrs are estimated. Hence, a very crude range for 
the precession period of these remarkable flows is 2800- 
14,600 yrs. 



3.4 Photometric characteristics of outflow sources 

Fig- E] demonstrates how the outflow sources differ from the 
bulk of the c2d-catalogue YSOs in Perseus. The outflow 
sources, plotted with black filled bars, are generally much 
redder and have a steeper near-IR/mid-IR SED than the 
general YSO population. The flow sources identified in this 
paper have a mean spectral index (a) of 1.39±1.14; the c2d 
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Figure 8. H2 2.122 fim image of the region to the south of NGC1333. Labelling and symbols are as in Fig|3] 



YSOs, which are predominantly Class II T Tauri stars, have 
a mean a of -0.32±0.89. (In the c2d catalogue, in regions 
with A v < 5 where sources are mostly background main 
sequence stars, the average value of a is -2.8: Evans et al. 
2005.) 



The outflow sources also occupy a relatively distinct 
region in colour-colour space (Fig. I10p . They possess a 
high [4.5]- [24] IRAC/MIPS colour, as predicted for young 
sources with ma ssive accretion disks IjWhitnev et al.ll2004l ; 
I Allen et al.ll2004h . However, they also exhibit positive [3.6]- 
[4.5] colours, unlike the bulk of the c2d YSOs, where the 
near-zero [3.6]- [4.5] colours are due to the intrinsic shape 
of the YSO SED. In Fig. [TO] 80% of the outflow sources 
that have measured 3.6 /im and 4.5 /im magnitudes have 
a [3. 6] -[4. 5] colour greater than 1. This is compared to only 
11% of the total c2d catalogue of YSOs in Perseus. Of the 
YSOs with a [3.6]-[4.5] colour greater than 1.0, 18% have 
outflows; only 0.4% of YSOs with [3.6]-[4.5]<1.0 drive out- 



flows. All factors point to the association of outflows traced 
in H2 with the youngest Class 0/1 protostars. 

In Fig.JTJJwe plot a against [3.6]-[4.5] colour for the en- 
tire Perseus c2d catalogue. Sources with values of a derived 
from just the two Spitzer bands plotted on the x-axis (i.e. 
those detected in only bands 1 and 2) lie along a narrow, 
diagonal strip. Sources detected in multiple Spitzer and/or 
2MASS bands exhibit some scatter about this strip, since a 
is derived from photometry points that are modified by var- 
ious absorption features (ice, silicate dust, etc.) and emis- 
sion lines (H2, [Fell], [Nell], etc.), as well as the intrinsic 
shape of the YSO SED (i.e. that of an embedded source 
with a disk and probably also a disk-hole). The Spitzer 
YSOs tend towards higher values of a and redder [3.6]- [4.5] 
colours, though they are largely indistinguishable from the 
other sources (field stars, background galaxies, etc.) in this 
plot. However, the outflow sources stand out as having ex- 
cess emission at 4.5 /im and thus red [3.6]- [4.5] colours. This 
is partly due to the inclusion of emission lines associated 
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Figure 9. Histograms of IRAC [3.6]-[4.5] colour and source spectral index, a, for YSOs in Perseus. Open bars represent all c2d identified 
YSOs; filled bars are sources with identified outflows. 




Figure 10. Spitzer colour-colour diagrams for c2d YSOs in Perseus (open circles). Filled circles represent sources with identified outflows. 



with infall and outflow in their integrated 4.5 fim fluxes. The 
4.5 /im Spitzer band covers the H2 pure-rotational 0-0S(9), 
S(10) and S(ll) lines at 4.69, 4.41 and 4.18 /im, as well as 
the CO v=l-0 band between 4.45-4.95 /im and the HI Bra 
line at 4.052 /tm, both of which may be excited in accretion 
flows. 

Some of the sources identified in the c2d catalogue as 
YSOs, based on their Spitzer colours and/or spectral in- 
dex, will in fact be molecula r hydrogen emis s ion-lin e knots. 
Shock models developed by ISmith fc Rosen! (|2005l ) predict 
[3.6]-[4.5] colours in excess of ~0.5, as well as positive val- 
ues for a (when derived from all four IRAC bands). These 
predictions are supported by recent Spitzer observations of 
HH 46/47, where the j et and bow shocks are detected in 
all four IRAC channels l|Velusamv. Langer fc Marshll2007l ) - 
note that there are pure-rotation H2 emission lines in all of 
the IRAC bands. 

In Fig.[TT]we plot the [3.6]- [4.5] colours and correspond- 
ing spectral indices for six shock models. The arrows indicate 



the effect of 5 magnitudes of K-band extinction (reddening 
moves the points to the upper-right in the figure). The ar- 
row labelled "Planar" represents a slab of molecular gas at a 
temperature of 2000 K with a density of 10 5 cm -3 . Arrows 
labelled C80, C50, C40 and C20 are C-type bow shocks: all 
have the same shape and post-shock density (10 s cm -3 ), dif- 
fering only in shock velocity (80 km s _1 , 50 km s , 40 km 
s" 1 and 20km s" 1 respectively). Decreasing shock veloc- 
ity pushes the C-bow data points toward the upper-right in 
Fig- EH Higher values of a and redder [3.6]- [4.5] colours are 
produced by a relative increase in emission from the lower- 
energy pure-rotational H2 lines (0-0S(4) and S(5) in IRAC 
band 4 and 0-0S(6) and S(7) in band 3) over the higher en- 
ergy 0-0 lines (S(9)-S(ll) in band 2 and S(13) and above in 
band 1) and ro-vibrational lines (1-00(5) and 0(7) in band 
1) at shorter wavelengths. The fastest C-type bow, C80, has 
a dissociative cap and is much like the J-type bow shock. It 
is therefore not too surprising that it coincides almost pre- 



10 C.J. Davis et al. 



m 3 













» 














/ 






/- 

- 










/ — 














C4C 


t 














C50 . 

C80 \ 




• * 










150 v 


/ ■ 




» • * 




3 










IP' 1 


Planar 




- 


















- 


•e _ 














japs a 


» 










* 
















* 


?3P 


1 











-2-10 1 2 3 

[3.6HL5] 



Figure 11. Plot of spectral index (a) against [3.6]-[4.5] colour. 
Grey dots represent the entire c2d Perseus catalogue, red open 
circles represent c2d identified YSOs, and black filled circles rep- 
resent outflow sources. Green circles mark c2d YSOs that may ac- 
tually be misidentified H2 knots. The arrows show the predicted 
positions for H2 emission from a Planar shock, a J-ty pe bow shock 
and four C-type bow shocks IjSmith &: Rosenll2005l ). The arrows 
reflect an increase in K-band extinction of 5 magnitudes. 



cisely with the 50 km s 1 J-type bow model arrow, J50, in 

Fig. ED 

The shock model data points overlap quite significantly 
with the outflow sources in Fig. [TT] The c2d "YSOs" that 
are most likely to be H2 knots, features 7, 12, 13, 16, 38, 
49 and 72 (the green circles in Fig. \\\\ have relatively low 
values of a and therefore appear to be associated with either 
fast C-type or J-type bow shocks. The low spectral indices 
measured for the observed knots could be caused by CO 
bandhead emission in IRAC band 2 or the inclusion of 24 /im 
data when calculating a for the observed knots. We do not 
include this data point in the model points, since there are 
no H2 lines in the MIPS 24 /im bandpass. However, weak for- 
bidden emission from [Fe ll] at 24.51 /im and 25.98 (like 
that seen in HH 46/47; IVelusamv. Langer fc Marshl 12007m 
could result in the "weak" detection of an outflow feature 
in this band in the c2d catalogue, and hence a lower value 
of a based on a fit to IRAC and MIPS data. CO emission is 
unlikely to contribute to observed outflow photometry data 
since very high densities, in excess of 10 6 cm -3 , are required 
to produce appreciable flux from this molecule. However, 
under such circumstances (seen perhaps in high-mass star 
forming regions) CO could account for as much as 10% of 
the observed emission in IRAC band 2. 



3.5 Outflows and the ambient medium 

The number of outflows identified here matches (to within 
20 %) the number of prot ostellar cores in L1448, L1455 and 
Bl jHatchell et alj|2007bj ). though it is lower by a factor of 2 
in NGC1333. In this densely-populated region this discrep- 
ancy is partly due to H2 features being identified as part of 
a single flow, when in fact they are associated with multi- 



ple flows. Even so, the association of protostellar cores with 
molecular (H2/CO) outflows is clear. 

It is also worth noting that the more massive clouds 
produce more outflows in Perseus- West. In Bl, NGC1333, 
L1448 and L1455 we identify (at least) 8, 11, 4 and 5 
outflows, respectively. From near-IR extinction mapping 
of Perseus , maps which cover an A v range of about 3-7, 
iKirk et all (|2006h identify large scale structures which they 
refer to as "supercores" . They stress that the bulk of the 
cloud material is contained in these extensive, moderate- 
density regions, rather than in the high-extinction cores 
not traced by the extinction map. Kirk et al. measure 
masses of 441 M , 973 M , 174 M Q and 240 M Q for the Bl, 
NGC1333, L1448 and L1455 supercores, respectively. Thus, 
in Perseus- West there is one molecular outflow for every 44- 
88 Mq of ambient material. Roughly 50 generations of ac- 
cretion/outflow activity would be needed to convert all of 
this molecular material into solar-mass stars. 

In some regions many of the jets appear to be roughly 
parallel, though this common angle is not obviously related 
to the large-scale cloud structure. L1448 is perhaps the best 
example, where the position angles of the four outflows iden - 
tified in Fig.|6]are between 147° and 155°. lKirk et~ail (|2006h 
fit 2-dimensional Gaussians to their supercores and, from 
their Fig. 5, we measure a supercore position angle of 120° 
for L1448 . In L1448, the outflows thus appear to flow along 
the length of the supercore. However, the two parallel flows 
in the Bl-ridge region (Fig.0] PAs of 124° and 130°) are al- 
most perpendicular to the elongated Bl-ridge supercore (PA 
~ 50°). The flows in Bl, NGC1333 and L1455 are more ran- 
domly orientated, although at least four of the larger flows 
in NGC1333 are aligned roughly north-south, parallel with 
the associated supercore (PA ~ 10°). Even so, we conclude 
that in general the outflows are randomly orientated and 
that the large-sc ale could structur e, as traced by the extinc- 
tion mapping of lKirket ail (T2006). has no obvious influence 
on flow orientations. 



3.6 Momentum injection and cloud support 

Molecular outflows from young stars inject turbulent energy 
and momentum into the surrounding medium. Over time 
these flows probably provide turbulent support against cloud 
collapse and may even hinder the star formation process. 
In low-mass star forming regions usually only a dozen or 
so outflows are observed at any one time (NGC1333 is a 
good example). Here we investigate whether these outflows 
provide sufficient kinetic energy to account for the observed 

turbulence in the ambient medium. 

As was noted by IWalawender. Bally fc Reipurthl 
(2005a), one can estimate the turbulent momentum and en- 
ergy in a star- forming region from the cloud mass, M c , and 
velocity dispersion, SVtuv, measured from (sub)millimetre 
molecular lines: 

• For Perseus - West we again use cloud masses from 
IKirk et all (2006). The masses of their supercores (estimated 
from 2MASS extinction mapping for Ay >3) are listed in 
Sect. 4.4. Although a sizable fraction of the overall cloud 
mass is likely to be in diffuse, lower- Av regions, outflows 
are unlikely to impart much of their momentum into these 
extended regions. Also, masses derived from an extinction 
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map are preferred over masses estimated from molecular 
line observations, since the latter may suffer from depletion 
and/or opacity effects. And like submillimetre dust contin- 
uum maps, optically thin molecular lines will only trace the 
denser core regions (Av >5). 

• A turbulent velocity is difficult to generalise for each 
region. The widths of molecular lines will result from tur- 
bulent and thermal motions in the molecular gas, as well as 
macroscopic motions if many independent cores are viewed 
along the same line of sight. In nearby, low-mass star form- 
ing regions the latter should not be a major problem. How- 
ever, due to opacity and depletion effects, different molecu- 
lar lines will trace different comp onents of the ambient gas. 
iKirk. Johnstone fe Tafailal < |2007f > find that in Perseus N 2 H+ 
(1-0) lines are broadened almost exclusively by thermal mo- 
tions in dense cores; C ls O (2-1) is more sensitive to lower 
density material in the envelopes of dense cores (the C 18 
transition has a critical density of 10 3 cm -3 and freezes out 
onto dust grains at densities above 10 cm" ~ 3 ). C 18 Ois there- 
fore a better tracer of the non-thermal motions in the inter- 
stellar - and intercore - environment pervaded by molecular 
outflows. Across Perseus C ls O profile widths vary by only a 
few tenths of a km s _1 . From Kirk et al. (2007) we there- 
fore adopt a mean value of 0.6 km s _1 for SVtm for all regions 
considered. 

In Table [5] we list the turbulent momentum P tur ~ 
McSVtui and turbulent energy E tU i ~ M c 5V t 2 UI in each of 
the clouds observed in Perseus. 

The momentum and energy of an outflow from a young 
star is hard to measure from infrared observations. There- 
fore, to estimate the combined momentum and energy, Pjets 
and Ej e ts, of the outflows from the embedded YSOs in each 
supercore in Perseus, we use the canonical values for the mo- 
mentum flux for Class and Class I low mas s protostellar 
outflows measured by iBontemps et al.l (|l996T ). and multi- 
ply these by the accretion lifetime of the embedded phas e 
jWalawender. Bally fc Reipurthl l2005al ; iDavis et al.l 12007ft . 
The momentum and energy flux of an outflow is thought to 
decr ease with tim e, in line with decreasing accretion rates 
(e.g. ISmithll2000h . We therefore ignore the Class II/III out- 
flow phase. If the momentum flux for the Class and I 
phases is ~ 6xl0~ 5 M© km s _1 yr _1 and ~ 4x 10~ 6 M km 
s _1 yr _1 respectively, then for durations of 10 4 and 10 5 
years, over its lifetime an outflow from a low-mass proto- 
star should inject ~1.0 Mq km s _1 of momentum into its 
surroundings. The outflow kinetic energy scales with the ve- 
locity; since most of the molecular material in outflows is 
probably entrained in the relatively slow-moving wings of 
bow shoc ks, then for a molecul a r flow velocity of ^20 km 
s- 1 (e .g. iKnee fc Sandell US iHatchell. Fuller fc Richerl 
l2007aft . the kinetic energy supplied by a typical flow would 
be ~ 4 x 10 37 J . In Table [2] we thus multiply these canonical 
values by the number of observed Eb flows in each region to 
get Pjets and Ej eis . 

Although the combined energy of the outflows, £j e ts, 
matches our estimates for the turbulent energy in each cloud, 
Etui, the combined momentum, Pjets, falls short of P tur by 
at least an order of magnitude in each region in Table [2] 
It has been argued that a large fraction of the jet mo- 
mentum may be carried by a collima ted, high-velocity jet 
|Walawender. Bally fc Reipurt h 2005a). However, the trans- 



port of this momentum component to the surrounding cloud 
may be relatively inefficient, as the heavy jet punches its 
way through the ambient medium. The existence of parsec- 
scale jets and the high proper motions of even the most dis- 
tant jet knots suggest that jets are relatively unhindered by 
their surroundings. Our assumption, that only the massive 
though low-velocity molecular flow component contributes 
to the turbulent motions in the ambient medium, is there- 
fore probably still valid. 

Instead, the turbulent momentum in the surrounding 
cloud may be built up via many generations of outflows. The 
lifetime of a Giant Molecular Cloud is 1-2 orders of magni- 
tude longer than the duration of the Class 0/1 phase, and 
Class 0, I and Class I I YSOs clearly co-exist in many star- 
form i ng regions (e.g . Gutcrmuth ct al. 2004; Kumar ~et al.l 
120071 ; lHarvev et all 120071 ; iGutermuth et al.1 l200Sft . Hence, 
over the lifetime of a GMC, many generations of outflows 
may contribute to the observed turbulent momentum in the 
surrounding cloud. 



4 CONCLUSIONS 

We examine near- and mid-IR images of a ~4 square de- 
gree region in Perseus- West and compare these to published 
catalogues of HH objects and moderately-extensive CO 3- 
2 maps. We find that in most cases molecular shock fea- 
tures are morphologically the same at 2.12 /j,m and 4.5 /im, 
although the more embedded portions of some flows are re- 
vealed only in the longer-wavelength data. H2 features are 
closely associated with high-velocity CO lobes, in support 
of bow shock entrainment scenarios. However, they often do 
not coincide with their optical HH counterparts, because of 
extinction effects or the discordant excitation requirements 
of each tracer. 

We use IRAC and MIPS photometry to identify and 
characterise the outflow sources. In comparison to the full 
sample of YSOs in the c2d catalogue, the molecular outflow 
sources possess extreme values of near-to-mid-IR spectral in- 
dex, q, and highly reddened [3.6]-[4.5] and [4. 5] -[24] colours, 
as is befitting their extreme youth. 

We find no correlation between H2 flow length and a, 
nor between flow opening angle (as traced in H2) and a. 
In general, the outflows are randomly orientated. However, 
we do find a correlation between the number of outflows 
and the number of protostellar cores in three of the four 
regions studies; in terms of total cloud mass, there is one 
molecular outflow for every 44-88 Mq of ambient material. 
The outflows may also be an important source of turbulent 
energy in the interstellar medium. 
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Table 1. Embedded YSOs and their molecular outflows in Perseus- West 
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"Embedded YSOs in Perseus- West; list compiled bv ljorgensen et alj l|2007h from Spitzer and SCUBA observations (see their Table 3). 

b H2 features, molecular outflows, and HH objects associated with each YSO. See Table [All for additional details. 

c Source photometry obtained directly from the c2d catalogues. 

d Mean flow length, position angle and opening angle (average of two lobes). 

e Unresolved YSOs associated with three separate, clearly defined flows (see the online Appendix for details). 
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Table 2. Turbulent momentum and energy in Perseus: clouds 
versus outflows 
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4 


1.6 


L1455 


144 


1.7 


5 


5 


2.0 



"Turbulent momentum and energy associated with the molec- 
ular cloud. 

6 Number of H2 flows in the region. 

c Total momentum and energy associated with the "N" 
molecular outflows in each region. 
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APPENDIX A: DISCUSSION OF OUTFLOW 
ACTIVITY IN PERSEUS- WEST 

In this section we present colour-composite images of each 
region constructed from WFCAM and IRAC observations, 
with 2.12 pm, 3.6 pm and 4.5 pm data occupying the blue, 
green and red channels, respectively (Figs IA1HA5|) . These 
show how UKIRT /WFCAM and Spitzer/IRAC data com- 
plement each other in studies of outflows. H2 features ap- 
pear pink/purple, since line-emission is in almost all cases 
detected at 2.12 /im and 4.5 pm. Higher-excitation H2 knots 
will appear slightly bluer, since the WFCAM filter covers a 
rotation- vibrational line while the IRAC band 2 covers only 
pure-rotational transitions (discussed further in Sect. 4.3). 
The outflow sources, many of which are not visible in the 
near-IR, are in most cases detected in the IRAC bands. The 
most embedded sources appear red in the colour images, as 
do the most deeply embedded outflow features that are not 
detected at 2.12 pm. 

These WFCAM/IRAC colour images are preferable to 
figures created from three of the four IRAC bands since one 
avoids using the less sensitive and relatively poor resolu- 
tion IRAC bands 3 and 4 data: in colour composites made 
solely from IRAC data (e.g. 3.6/4.5/8.0 pm) faint, embed- 
ded sources may be undetected at the longer wavelength and 
thus will appear green, similar to H2 line-emission features 
that are bright in IRAC band 2. The higher resolution WF- 
CAM H2 2.122 pm data also help to "sharpen" the outflow 
features, producing an effect similar to unsharp-masking. 

Below we discuss the outflows in each region separately, 
relating the H2 features detected at 2.12(i.m and/or 4.5(J.m 
with the high- velocity CO outflow lobes displayed in Figs. [3}- 
[7] the HH objects listed in the literatu re, and the proto- 
stars identified from S pitzer photometry |Evans et al.ll2005l ; 
Ijprgensen et alj|2007h. Note that e mbed ded YSOs from the 
list presented by J0rgensen et al.l (|2007T l are referred to as 
YSO 1, YSO 2, etc., while sources from the c2d catalogue 
are identified by their RA and Dec, and are prefixed by 
"J ". Alternative s ource names are obtained from Table Al 
of lHatchell et all (|2005f ) and refere nces therein; NGC133 3- 
IRAS source names are taken from I Jennings et al.l (|l987T ). 



Al Bl 

We identify six outflows in Bl. All six are driven by sources 
that are ide ntified by c2d and all ex cept LkHa327 are in the 
catalogue of |j0rgensen et all l|2007h . (There is potentially a 
seventh outflow that may have a source that is undetected 
by Spitzer between Bl and HH 429 - see below.) 

The deeply embedded Class source YSO 38 (SMM 1, 
Bl-c) is associated with a well-defined bipolar CO outflow 
and the H2 features 5 and 6-8 (Fig. [3]). The H2 data show 
that the outflow extends well beyond the bounds of the 
CO map. In our colour WFCAM/Spitzer image (Fig. IA1|) . 
an 'S-shaped', precessing flow is evident; the flow is traced 
back to the source in the longer wavelength 4.5 pm emis- 
sion. YSO 38 clearl y drives knot 5 and prob ably much of 
object 6-8 (see also Wal awender et ai"1l2005bT ). Parts of the 
6-8 complex could be powered by J033309. 3+311012 (the 
YSO situated midway between features 6 and 7 in Fig. [3]), 
though this source has a flat spectral index (~0.08) and 
only a moderately red [3.6] - [4.5] colour (^0.80; see e.g. 



Fig. [2}. This source may instead be a bright outflow feature 
misidentified by c2d as a YSO candidate (discussed further 
in Sect. 4.3). 

The overall morphology of the H2 complex 6-8 also adds 
to the confusion. Its east-facing bow shape seems to sug- 
gest movement toward the east. 6-8 could conceivably be 
a bow shock in the counter-lobe of a flow that drives the 
bright, opposite-facing HH object HH 429 (H2 knot 53; see 
Walawender et al. 2005b - Fig. 30). HH 429 is 17.5' due west 
of feature 53. 

To the south of YSO 38, the nebulous, Class I source 
YSO 37 (Figs. |3] and IA1|) powers knot 10 to the west and a 
red-shifted CO lobe to the east. This east- west orientation 
supports the inclusion of HH 787, 785, 781 (and possibly 
HH 429) to the west and HH 789 to the east in another 
parsec-scale flow (| Walawender et a l. 2005b). 

The blue-shifted CO emission to the southwest of 
YSO 37 is probably associated with a bipolar flow from 
YSO 35 (SMM3), although no H2 emission is detected from 
this flow at 2.12 pm. 

YSO 39 (Bl-b , SMM2,3) compris es two Class sources 
(Bl-bN and Bl-bS; lHirano et al.ll999l ). YSO 39 is associated 
with a bright, knotty H2 flow 9 to the southwest, and a blue- 
shifted CO lobe to the northeast. 

YSO 40 is a nebulous infrared source tentatively as- 
sociated with blue-shifted CO emission and faint H2 line 
emission (feature 50) . The source has a relatively steep spec- 
tral index, as expected for an embedded outflow source (a 
= 2.54, [3.6]-[4.5] = 2.02). It is possibly associated with 
the bright HH features H H 431,433 and the jet HH 790 
l|Walawender et al.ll2005bl , see their Fig. 32). 

Finally, LkH«327, the brightest near-IR source in Bl, 
drives a northwest -southeast outflow and ex cites HH 432 
and H2 feature 63. IWalawender et all (2005b) suggest that 
HH 430 and HH 788 trace the extended regions of the north- 
ern lobe, with HH 432 and 63 (their MH 7) tracing the 
southern lobe. 



A2 Bl Ridge 

The ridge of molecular gas running from Bl to L1455 
contains two J0rgensen YSOs that are clearly driving 
molecular outflows (Figs. gj. YSO 33 (SMM 76, IRAS 
03292+3039) is associated with H2 features 2, 3 and 4, 
while YSO 32 (SMM 77, IRAS 03282+3035) powers H 2 fea- 
tures 49 (HH 773 ) and 50. Bow shock 49, first noted by 
iBallv et al.l (I1993T ). is bright in H 2 2.122 pm (Fig. but 
only marginally detected in the optical, as expected for em- 
bedded flows from Class objects. The YSO 33 flow has no 
known HH counterparts. 

51 and 52 are H2 features located southwest along the 
Bl ridge (see Table EH for details). The area has several 
Spitzer-identified YSOs and red sources, along with several 
HH objects which include HH 368-372, 427, 428, 769, and 
792. However, none o f the YSOs are in the catalogue of 
|j0rgensen et all (120071) (even thoug h SCUBA data do ex- 
ist for this region: Kirk et aLll2006T ). The HH objects and 
H2 features do not have any clearly identifiable sources, al- 
though H 2 features 51, 52 and IRAS 03272+3013 do lie on 
the same axis. 
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Figure Al. SEE GIF UKIRT H2 2.122 fim (blue) - Spitzer 3.6 Aim (green) - Spitzer 4.5 (im (red) colour composite image of Bl. H2 
features are numbered; embedded YSOs from ljorgensen et alj i2007h are labelled. 

Figure A2. SEE GIF UKIRT H2 2.122 ^tm (blue) - Sp itzer 3.6 (green) - S pitzer 4.5 /im (red) colour composite image of Bl Ridge 
(west). H2 features are numbered; embedded YSOs from j0rgcnscn ct al. (2007) are labelled. 



A3 L1455 and L1451 

In the L1451/L1455 region in Fig.[5]we identify five H2 out- 
flows. At least four of these are driven by embedded YSOs 
l|j0rgensen et aL 1 12007!) . 

The most prominent outflow emerging from the L1455 
core includes HH 279, H2 knots 11-13 and possibly HH 423 
and H2 feature 85 in a parsec-scale flow (see Table [TJ. Class 
I objects YSO 5 (SMM 39) and YSO 8 (RNO 15) are can- 
didate sources, although neither is associated with high- 
velocity CO (Fig. [5). The c2d source J032719.0+301718 (co- 
incident with HH 279) has a red [3.6]-[4.5] colour of 1.55, 
though it also possess a negative spectral index (-0.46) and 
therefore could again be a H2 knot mis- identified as a YSO. 
YSO 8 may instead drive a north-south flow that includes 
H2 knots 14, 54 and weak blue-shifted CO emission to the 
south. 

Between YSOs 5 and 8 in the southeast corner of Fig. [S] 
YSO 6 (SMM 35, RN015-FIR, IRAS 03245+3002) powers 
H2 features 15-17. 55 and a precessing bipolar CO outflow 
<|Davis et al.lll997f ). YSO 7 (SMM 36) drives a second, or- 
thogonal CO flow that includes the faint H2 feature 59. The 
blue lobe of this CO outflow may extend ~1' northwest of 
knots 15 and 16, although no H2 features are detected in this 
extension. In the colour-composite image in Fig. IA3I YSOs 
6 and 7 are obvious as the reddest sources in the region. 

To the west of L1455, toward L1451, there is a parsec 
scale flow driven by YSO 4 (SMM 80). The flow is associ- 
ated with H 2 features 56, 57 and 58 (HH 492, 493 and 317, 
respectively) to the west, and the l ong chain of optical knot s 
that comprise HH 280 to the east (W alawender et alj [2004). 
Red-shifted CO extends towards this eastern shock feature. 



A4 L1448 

There are four outflows identified in Fig. [6] The reddened 
source YSO 1 (SMM 30, IRAS 03222+3034, L1448-IRS2) 
drives a spectacular bipolar flow associated with H2 fea- 
tures 72 and 74 (Fig. El see also Eisloffel 2000). Additional 
bright line emission features are seen to the northwest of 
YSO 1, although these are probably associated with a sec- 
ond bipolar flow from L1448-IRS1 (RNO 13). L1448-IRS1 
is not catalogued as a YSO by j0rgensen. It has a spectral 
index of 0.38. It is detected only in the Spitzer MIPS 1 band 
in the c2d catalogue and is classified as "star+dust". The 
lack of a detection in the IRAC bands is probably due to 
saturation as IRS1 is quite bright (see Fig. I A4p . H2 knot 81 
and HH 194 occupy the counter-lobe of this flow. 

The parallel outflows from YSO 2 and YSO 3 
(SMM 28/L1448N and SMM 29/L1448C, res pectively) have 
been studied extensively i n the past (e.g iBachiller et al.l 
Il99d ; iDavis fc Smith! 1 19951 : lEisloffell 120001 and references 
therein). The flow from YSO 2 includes H2 features 82 to 
the southeast and 78, 79, 80, 83 and 84 to the northwest. 



These features increase the length of this northwestern flow 
lobe to 25.5' (2.2 pc). The YSO 3 flow is associated with 
H2 object 77 in the south and the spectacular, curving H2 
flow 76 to the north. Feature 75 is an unusual group of knots 
to the northwest of 76: we tentatively associate these with 
the curving 76 outflow, though on morphological grounds 
one might regard 75 and 78 as bow shocks and counter-bow 
shocks in a bipolar flow centred midway between these two 
objects. However, no candidate source is found in this re- 
gion. 

Finally, we bring the reader's attention to a highly red- 
dened source 32" south of YSO 3. This nebulous source can 
be seen clearly in the colour image in Fig. IA4I The source, 
J032539. 6+304336, is classified by c2d as "cup-down", i.e. 
it has been detected in only 3 IRAC bands with a magni- 
tude at 4.5 /im that is brighter than at 3.6 /im and 5.8 ^tm 
(possibly due to unresolved line emission). However, close 
examination of the Spitzer IRAC images reveals that the 
source is in fact detected in all four IRAC bands. From our 
own photometry, measured from IRAC BCD images using 
a 2" aperture, we measure a [3.6]- [4.5] colour of 2.2 (typical 
for outflow sources) and a spectral index, a, of 2.33 (again 
typical for embedded Class 0/1 outflow sources). Hence, this 
source could be contributing to the outflow activity in this 
region. 

A5 NGC1333 

The outflows in NGC1333 have been studied in considerable 
detail, though the WFCAM images presented here are the 
most extensive H2 images to date and as a result reveal new 
features. In Fig. [7]and[8]we ma rk 11 outflows. Some of these 
have previously been noted by Hodapp fc La dd (1995]). All 
11 are associated with I j0rgensen et al. ( 20071 ) YSO sources. 



A 5.1 The SVS-13 SMM core 

There are three major outflows emanating from the gen- 
eral area of SVS 13 a nd YSO 20. YSO 20 is identified by 
|j0rgensen et al.l l|2007h as the centre of the SCUBA submil- 
limeter core rather than a YSO. This is because YSO 20 was 
identified using only SCUBA data and so no c2d YSO po- 
sition is available. This s ubmillimeter core conta ins 5 radio 
continuum sources which iRodriguez et al.l l| 19991 ) identify as 
YSOs. 

Rather than speculate on the source of each of the flows 
emerging from the SVS 13 core, in Table [TJ we identify the 
sources according to their alignment with outflow features 
and/or CO flow lobes: YSO 20a is deemed to be the driver 
of the well-know bipolar CO outflow containing HH 7-11; 
YSO 20b is the driver of the north-south flow at PA ~ 170° 
which includes H2 features 30 in the south, 37 in the north 
and possible feature 38 (HH 12); YSO 20c is then the driver 
of the lone south lobe at PA ~ 10° . This latter chain of H2 
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Figure A3. SEE GIF UKIRT H2 2.122 /im (blue) - Spitzer 3 . 6 fim (green) - Spitzer 4.5 /im (red) colour composite image of L1455. H2 
features are numbered; embedded YSOs from j0rgen sen et al. I (|2003) are labelled. 



Figure A4. SEE GIF UKIRT H 2 2.122 /im (blue) - Spitzer 3 . 6 ^m (green) - Spitzer 4.5 (im (red) colour composite image of L1448. H 2 
features are numbered; embedded YSOs from |j0rgensen et al.l J2007I) are labelled. 



features extends at least 16.5' (1.4 pc) to the south and in- 
cludes H 2 /HH objects 34, 28, 43, HH 14 and 19. There is no 
obvious counter-lobe in our H2 image, although HH 754/755 
do lie north along an axis drawn through the southern lobe 
i|Walawender et al. 2005b). There is also no clear evidence 
of an outflow in the CO data of lKnee fc Sandelll (2000). Our 
colour-composite image in Fig. IA5l reveals a bright, nebulous 
source roughly 1' south of SVS 13 that could be associated 
with the YSO 20c flow. The source is associated with a con- 
ical reflection nebula that opens out toward the south. 

A5.2 IRAS 2 A and 2B 

YSOs 15 and 16 are better known as IRAS 2A (SMM 44) 
and IRAS 2B. YSO 1 5 is associated with a SCUBA core 
and VLA radio source (jRodriguez et al] 1999). In the colour 
image in Fig. IA5I YSO 15 appears highly reddened and elon- 
gated along the identified flow direction. It is clearly asso- 
ciated with H2 knots 31 and 32, and the bipolar CO flow 
evident in Fig. [7] However, the YSO 15 "core" may also be 
associated with a second, orthogonal flow, which we label 45. 
This bipolar flow can also be seen clearly in the CO data. 
The north-south flow from YSO 15 is a large outflow that 
influences HH 12 in the north and extends at least as far 
down as HH 13 and H2 features 41/42 in the south. Its axis 
also lines up with knot 18 (HH 746) ~25' to the south of 
YSO 15. The extent of this flow lobe could therefore be at 
least 2.2 pc. 

There are no H2 features that are obviously associated 
with YSO 16 (IRAS 2B). 

A5.3 HH 12 

HH 12 is a complex cluster of knots about 5' north of SVS 13 
and HH 7-11. The bright infrared source (IRAS 6) associated 
with HH 12 sits atop a colum n of dense, molecular g as; the 
source is flagged as a YSO bv I Jeirgensen et all l|2007h based 
on its association with a SCUBA core (SMM 45). 

We label the H 2 features around HH 12 as 37/38. Fea- 
ture 37 is a faint trail of emission leading from near SVS 13 
to the head of HH 12; 38 is the complex group of bright knots 
that make up the head of HH 12. Part of HH 12 could be 
associated with the H 2 flow 31/32 from YSO 15/IRAS 2A; 
HH 12 is clearly aligned with the well-defined axis of this 
bipolar outflow. However, the spray of knots to the west of 
YSO 19, coupled with the overlapping blue-shifted and red- 
shifted CO observed towards this source, point to an outflow 
inclined towar ds the line of sigh t that is driven by a source 
within HH 12. iRodriguez et all (|l999h find 2 radio contin- 
uum sources that coincide with the YSO 19 core and there 
are three flat-spectrum Spitzer YSOs about 1' to the south- 
west of YSO 19 that could be contributing to the outflow 
activity, although one or more of them could again be an H2 
knot misidentified as a c2d YSO. 



A5.4 Other outflows in NGC1333 

Northeast of SVS 13, H2 feature 36 (HH 6) and its associated 
CO outflow are driven by either YSO 23 or YSO 24 (IRAS 7, 
SMM 46). There is also an extensive flow perpendicular to 
36 that is likewise driven by either YSO 23 or 24. This large 
flow runs from H2 knot 44/HH 5 in the south to feature 68 
in the north (see Fig. \M^. Feature 25/HH 18 and HH 761 
may also be part of this flow. In Tables [Tl and I All YSO 23 is 
identified as the source of the extensive northeast-southwest 
jet and YSO 24 as the source of HH 6, because YSO 24 is 
slightly better aligned with the central axis of the feature 
36/HH 6 cluster. However, the reverse cannot be ruled out. 

NGC1333 also contains two "S-shaped" outflows, driven 
by YSOs 11 and 14 (Figs[7]and[8]|. YSO 11 (SMM 49) drives 
H 2 features 20, 46 and possibly 71/HH 338 (8.1' to the north- 
west), and feature 21 and possibly even 24/HH 15 (7.3' 
to the southeast). YSO 14 drives H 2 knots 22/23/62 and 
HH 340/343 in another S-shaped flow (for a detailed discus- 
sion of this outflow see iHodapp et ail [2005). Both flows are 
remarkably similar to the bending jet from YSO 38 in Bl 
(Fig.lAlJ. 

3' southeast of SVS 13, YSO 22 (SMM 41) drives a well- 
defined bipolar CO outflow that includes H2 features 33 and 
29 (Fig. [TJ. To the southwest, this flow may be interacting 
with the flow from the SVS 13 core and to the northeast 
with both the HH 7-11 outflow and the southern lobe of the 
YSO 23 jet. 

Finally, near the bottom of Fig. YSO 18 (SMM 65) 
drives an east-west flow associated with H2 bows 26 to the 
east and 28 to the w est. Unfortunately the CO map of 
iKnee fc SandelH l|2000l ) does not extend sufficiently south- 
ward to contain YSO 18 and its outflow, or indeed any of 
the other southern flows labelled in Fig. [8] 
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Figure A5. SEE GIF UKIRT H 2 2.122 /nn (blu e) - Spitzer 3.6 Mm (gre en) - Spitzer 4.5 yum (red) colour composite image of NGC1333. 
H2 features are numbered; embedded YSOs from j0rgenscn ct al. (20o3) are labelled. 
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a Position of the brightest knot within the H2 fe ature or flow. 

& Most likely source. Numbers refer to Table 3 in |j0rgensen et al ] (|2007h . A dash means there is no obvious source candidate. 
c The distance from the source to the most distant knot. If no known source, entire length of all knots. 
^Position angle, measured east of north, of an axis drawn from the source to the feature. 

e Associated HH object (if any is known). 

Prev ious names of H2 features. AS R refers to sources fromlAspin. Sandell. fc Russelll jl994h : HL refers to sources from lHodapp fc Ladd 
lll995h ; MH refers to sources from [Walawcnd er et al. I ll2005bf i 
9 Figures in which the feature appears. 
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